function result = hellinger_norm(p_vec, q_vec)
%HELLINGER_NORM Calculates Hellinger distance for discrete disterbutions -
%               p_vec and q_vec are the distributions.
    result = norm(sqrt(p_vec) - sqrt(q_vec),2)/sqrt(2);
end

